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The Casimir pressure is calculated between parallel metal plates, 
containing the materials Au, Cu, or Al. Our motivation for making 
this calculation is the need of comparing theoretical predictions, based 
on the Lifshitz formula, with experiments that are becoming gradu- 
£h ' ally more accurate. In particular, the finite temperature correction is 

considered, in view of the recent discussion in the literature on this 
q-i point. A special attention is given to the case where the difference 

between the Casimir pressures at two different temperatures, T = 300 
K and T = 350 K, is involved. This seems to be a case that will be 
experimentally attainable in the near future, and it will be a critical 
d ' test of the temperature correction. 
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Figure 1: Parallel plate-geometry. The two surfaces are located at z = and 
z = a. 

1 Introduction 

The Casimir effect pQ has in recent years attracted a great deal of interest 
(for recent reviews, see 13 HI EI)- The advent of accurate experiments 
has accentuated the need of performing detailed calculations of the Casimir 
forces, based upon realistic input values for the permittivities in the (assumed 
homogeneous) materials. In the case of two semi-infinite media separated 
by a gap a - the standard set-up in the Casimir context - the formula in 
question is that due to Lifshitz [6 . In the case of a micrometer-sized sphere 
above a plane substrate - a case that is tractable via use of the proximity 
force approximation when the spherical surface is weakly curved [Zj - the 
experimental accuracy is claimed in the literature to be on the 1 % level. 
We shall not here give an overview of recent experiments; we will return to 
some examples below. The reader may instead consult recent reviews: a 
detailed exposition on the experiments up to 2001 is given in Bordag et al. 
[I] , a survey of the developments in the last four years is given by Milton [2] , 
section 3.6, and the works of the Purdue group is presented by Decca et al. 
jSJ. (A brief survey of the experiments is given also in the note [§].) 

We shall consider the simple set-up shown in figure H There are two 
metallic semi-infinite media of permittivities E\ and £3, with a dielectric 
medium of permittivity e 2 in between. For simplicity we assume that re- 
gion 2 is vacuum (air), so that e 2 = 1. The surfaces are assumed to be 
perfectly flat, of infinite extension, and the media are assumed nonmagnetic. 
Our intention is to work out values for the attractive Casimir surface pres- 
sure jF(a) versus gap width a for similar and dissimilar metals, when the 



temperature T is finite. Of main interest will be the temperature correction, 
in view of the conflicting opinions in the literature on this point. We will 
follow the same calculational strategy as in our earlier recent papers on these 

issues dunning. 

We shall consider three different metals; gold, copper, and aluminium. 
For these metals we have access to excellent numerical data for the permit- 
tivities (courtesy of Astrid Lambrecht and Serge Reynaud). We know how 
varies with imaginary frequency ( over seven decades, ( G [10 11 , 10 18 ] 
rad/s, at room temperature. For frequencies up to about 1.5 x 10 15 rad/s 
the data are nicely reproduced by the Drude dispersion relation 

where uj p is the plasma frequency and v the relaxation frequency. For the 
three metals mentioned we have |T3*1 ITI] 

u p = 9.0eV, z/ = 35meV Au 

u p = 9.05 eV, v = 30 meV Cu 

Wp = 11.5eV, z/ = 50meV Al (2) 

(note that 1 eV=1.519 x 10 15 rad/s). Using these data we can calculate the 
Casimir pressures to an accuracy better than 1%. 

We shall consider three different temperatures. First, it is of interest to 
work out explicitly the zero-temperature Casimir pressure. When discussing 
finite temperature corrections one should first of all know what is meant by 
the T = reference level. This issue is not trivial, since most of the T = 
theoretical predictions have been referring to the idealized case where e = oo 
from the outset. As discussed extensively in earlier works 0111111121111], the 
correct model in an idealized setting is the modified ideal metal (MIM) model, 
which assumes unit reflection coefficients for all but the transverse electric 
(TE) zero frequency mode. Our argument rests upon the condition that the 
relaxation frequency z/(0) at zero frequency remains different from zero. Here, 
as in Ref. ^1], we wm calculate the T « pressure numerically, inserting 
real data for e(i(). We shall choose T = 1 K as the lower temperature 
limit. It turns out that this limit is stable numerically, and numerical trials 
around this limit indicate that it describes the zero temperature case with 
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good accuracy. This method, although numerically demanding, is physically 
better than adopting the simple idealized metal model. 

The second temperature of interest is room temperature, T = 300 K. 
Recognizing that the difference between Casimir pressures at T = and T = 
300 K will hardly become a measurable quantity we shall instead consider, 
as our third chosen temperature, T = 350 K. The difference between the 
Casimir pressures at the two last-mentioned temperatures will perhaps soon 
become accessible in experiment. We shall therefore focus upon calculating 
how this difference varies with a. 

The present calculations, involving three different metals and three dif- 
ferent temperatures, are of interest for comparison of future experimental 
work on the finite temperature Casimir force. To our knowledge, these are 
the first calculated results for unequal metal surfaces at finite temperature. 
We ought here to mention that the first correct calculations of finite tem- 
perature Casimir forces between two equal gold, copper, or aluminium metal 
half planes were made by Bostrom and Sernelius JSJ ; these authors con- 
sidered also two thin equal metal films of gold, silver, copper, beryllium, or 
tungsten [T71[TB|. The calculations by Bostrom and Sernelius showed that 
the retarded van der Waals or Casimir interaction in more or less the same 
separation range (around 1 //m) as considered in the present paper depends 
on the choice of material for two equal metal surfaces. The present work is 
thus a step beyond this in that it considers unequal metal surfaces. 

It has repeatedly been pointed out by some authors that our use of the 
Drude dispersion relation runs into conflict with the Nernst theorem in ther- 
modynamics; cf., for instance, [T^HID!- We have shown earlier, however, that 
this is not the case fTJEI]. Thus, these thermodynamic issues will not be 
given further attention here. 

In the next section we present the general formalism, for similar as well as 
for dissimilar media, and give then in section 3 the results of our calculations 
in several diagrams. With respect to the temperature correction, we restrict 
ourselves in this section to the difference between T = and T = 300 K 
predictions. In section 4 we focus our attention on the difference between 
the mentioned 300 K and 350 K cases. 

In the main text, we put h — c — — 1. 
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2 Basic formalism 



We consider first the case of two identical media, E\ = e% = e. With the 
same notation as in ^UJ we can write the Casimir pressure as 



7r/3cr 
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(4) 

The prime on the summation sign means that the m = term is counted 
with half weight; (3 = 1/T is the inverse temperature. The minus sign in Eq. 
Q means that the force is attractive. 

[The following point should be noted. Assume that a plane wave is inci- 
dent from the left medium (z < 0) towards the boundary at z — 0. For the 
TM mode, the ratio between the reflected wave amplitude R™ and the in- 
cident wave amplitude A™ is equal to the square root of the coefficient A m , 
after the real frequency uj has been replaced with the imaginary frequency £ 
(uj = corresponding to the surface mode: 



R 



TM 



J^TM 



cf. Appendix B in ^Uj. Analogously, for the TE mode: 

R TE 
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(5) 
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Generalization to the case of dissimilar media leads to the following expres- 
sion: 
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where 



A TE _ S 1~P A TE _ S 3~P 
Sl+P s 3 +p 



A TM = £lP-Sl A TM = £ 3P~ S \ (g) 

eip + si - e 3 p + s 3 

Again, if the two media are equal we have Ai = A 2 for each of the modes, 
so that 

A^ M A™ — > A m , Aj E A^^B m , (9) 

and the formula (J3J) is recovered. 

We calculate the expression (JJJ) by means of MATLAB. The zero fre- 
quency case m — > may be treated separately by analytical methods, at 
least if we are considering an idealized model for the metal, because in this 
limit there is an interplay with the other limit e — > oo in the expressions for 
the coefficients. We have in this case p — > oo, — > p for % = 1,3, implying 
that 

A^^O, Aj M ^ £ -^ = A i . (10) 
Then, the m = contribution can be written as 



= -^/o, (11) 



where 



1 r 00 o , Ae- 2?/ I 



7 » = -2/o ^ r 3^ = -iP 'y'°s( 3 . A ). ( 12 > 

the polylog function being defined as polylog(x, z) = Y%?=i z n n~ x . We here 
assume that A = AiA 2 < 1. (For A > 1 the integral is undefined.) Since 
Si ^> 1 for a metal near ( = we have that Awl, but still less than 
unity, so we can let J — > — |polylog(3, 1) = — |C(3), where ((x) is the 
Riemann zeta function. The numerical value of Iq used in our calculations 
was I = -0.1502571129. 

The first equation in (jlOj) means that there is no contribution to the 
Casimir force from the m = TE mode, at finite temperatures. (At T = 
the m = effect vanishes, as the discrete Matsubara sum is replaced by 
an integral over frequencies.) This behaviour is a consequence of the Drude 
relation at low frequencies. The same behaviour can also be seen by use of 
quantum statistical methods, as was demonstrated for the case of spherical 
geometry in |2U E2] • 
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a = 50nm 




10° 10 1 10 2 

Dimensionless integration variable y 



Figure 2: Integrand of equation (J2J) versus y = qa for a = 50 nm, m = 
45000, T = 1 K. The materials are Al and Cu. Left end point of the curve 
corresponds to the lower limit y = vnq of the integral in ©. 



3 Numerical calculations. Results 

Since expression ((Zj) is complicated, with an upper limit y = qa = oo for 
the integral, it is useful first to get information about how the integrand 
varies with respect to y in typical cases. The expression is most demanding 
numerically for low temperatures and small gap widths. 

Figure El shows how the integrand varies with respect to y for a very 
high Matsubara number, m = 45000, when a = 50 nm and T = 1 K. The 
configuration is one aluminium and one copper plate. It is seen that when 
y becomes larger than about 10, the contribution to the integral decreases 
rapidly. At y = 50 the value is only 10~ 45 . Figure El shows for comparison 
how the integrand varies with y at the same temperature when the frequency 
is at the lowest non-vanishing value, m = 1, for a = 50 nm and a = 3000 
nm. The behaviour is seen to be quite similar to that above; for instance, 
when y = 50 the integrand becomes approximatively 10 -40 . For different 
plate combinations, we get approximately the same behaviour. For high 
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Figure 3: Same as figure but with a = {50, 3000} nm, and m — 1. 



temperatures, the same conclusion can be drawn. In all, we found it sufficient 
in our computations to adopt the value 

Umax = 50 + m7 (13) 

as general cutoff. [It is useful to note that 7 = 2744 (aT), when a is given in 
meters and T in degrees kelvin.] 

Numerically, we used a method of higher order recursive adaptive quadra- 
ture. This method approximates the value of the integral with a chosen 
tolerance of 10 -10 . 

Next, it is useful to show the variations of e(i() graphically for the three 
metals mentioned, together with information about the frequency region ac- 
tually used in the calculations. Figure |U shows this for the case of T = 300 K. 
The data extend over 7 decades. The vertical lines show that the important 
frequency region in this case lies between 2 x 10 14 rad/s and 5.5 x 10 16 rad/s. 
We also show the corresponding values of m = ( m f3/2ir (in nondimensional 
units). The range of m is [1,220]. For comparison, the theoretical predictions 
are shown for the case of gold, both when using the Drude relation (0) and 
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Figure 4: Numerical permittivity data for Al, Au and Cu (courtesy of Astrid 
Lambrecht and Serge Reynaud). Vertical lines show the frequency region 
that we used in the m summation at T = 300 K. Predictions from the Drude 
and the plasma dispersion relations are also shown in the case of Au. Top 
axis gives the values of the Matsubara number m. 



when using the plasma dispersion relation 

e(»C) = l + ^. (14) 

We see that for ( < 1.5 x 10 15 rad/s the Drude curve fits the data nicely, 
but for ( > 2 x 10 15 rad/s it gives too low values for e. The used frequency 
region corresponds to the area where the Drude prediction and the plasma 
prediction are approximately equal, and also to the area where the data for 
aluminium differ the most from the data for gold and copper. 

Figure 03 shows the analogous situation when T = 1 K. It becomes now 
necessary to use a much larger frequency region, from 7 x 10 11 rad/s to 
3.3 x 10 16 rad/s, corresponding to about 80% of the entire data set. The 
Matsubara number region is m G [1,40000]. The physical reason for this 
behaviour is, as always, that the case of low temperatures implies that the 
frequencies are very closely spaced. 
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Figure 5: Same room temperature data as in figure EJ but the vertical lines 
show the frequency region used in the calculations at T = 1 K. Drude and 
plasma dispersion relations shown for Au as before. 

3.1 Room-temperature Casimir force 

We shall in this subsection assume that T = 300 K. The force, expression (jJJ), 
is always negative, but it is convenient to represent it graphically in terms 
of the modulus |jF(a)|. Since in the parallel plate-experiment of Bressi et al. 

one was able to control the gap width down to 50 nm, we choose a = 50 
nm as our lower limit. As upper limit we choose a = 3 /im. To better visualize 
the results graphically, we divide all data sets into two groups, a e [50, 200] 
nm and a G [200, 3000] nm. As mentioned earlier, the plates are assumed 
infinite, and all roughness corrections are ignored. 
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Figures El and [7| show how the Casimir pressure varies with a, for various 
combinations of metal plates. The force between two gold plates was com- 
puted earlier in but is included here for comparison. The differences 
between the various combinations of the materials are seen to be small, and 
they diminish with increasing gap widths. The largest force always occurs 
for two aluminium plates. This may be called group I. The combinations 
aluminium-gold and aluminium-copper yield a somewhat smaller force (group 
II), and the last combinations Au-Au, Au-Cu and Cu-Cu result in the weak- 
est set (group III). When a increases from small to large values, the internal 
order in strength between the materials in groups II and III is interchanged. 

3.2 Room-temperature correction, compared to T = 

We now turn to the finite temperature correction in the Casimir force. As 
mentioned above, it is then important to make clear what we mean with 
the zero-temperature force. Numerically, we have seen that it becomes sat- 
isfactory to represent the latter case by the choice T = 1 K, the difference 
between T = and T = 1 K being negligible. 

Since the data for e(iQ are measured at room temperature, the natural 
question becomes: can we use these data also at very low temperatures? 
This issue has been discussed earlier, in [TTJ [T2J 123 , with the conclusion that 
the temperature dependence appears not to influence the dispersion relation 
in a way that changes the Casimir force significantly. This implies that we 
can insert the same Lambrecht-Reynaud data as before, and also assume the 
same values for u p and v (cf. equation (J2J)). 

The calculated results from the T = 1 K case are hardly distinguishable 
from the T = 300 K case when plotted in the same figure. Therefore, it 
is better to show the calculated finite-temperature differences. In figures |H1 
and 01 we show the difference between the Casimir pressures at T = 1 K and 
T = 300 K, 

AF = - IJ300I, (15) 

versus a in the range from 50 nm to 1700 nm. As the graphs for the various 
combinations of metals become indistinguishable for a > 0.4 fim and ap- 
proach zero when a becomes larger than 1.7 /im, we have omitted the region 
from 1.7 fim to 3 /im in figure 01 
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Figure 8: Difference between the Casjnair pressures at T = 1 K and T = 300 
K, equation (fTKjl . versus gap width a for a £ [50, 200] nm. 



An important property seen from the curves is that AjF is positive. That 
means, the force is weaker at room temperature than at T = 0. This is 
the same effect as we have pointed out earlier, in connection with identical 
materials in the plates [TJl EH HI] ] see also [21 [2B] • The behaviour is a direct 
consequence of the Lifshitz formula in combination with realistic permittivity 
data for the materials, the latter being, as we have seen, in agreement with 
the Drude relation at frequencies ( < 1.5 x 10 15 rad/s. Our results for 
the temperature dependence are in contrast to those obtained by use of the 
plasma dispersion relation; in that case, the deviation of the force is positive 
instead of negative, and is moreover very small [T^l 127]. 

To get an overview of the magnitudes of the temperature correction, let 
us give some examples: 

1. For small gap widths, the correction is relatively small. Thus when 
a = 100 nm, the Casimir pressure is 6.105 Pa at T = 1 K and 6.061 Pa at 
T = 300 K, thus giving a room temperature reduction of 0.72%. 

2. When a = 200 nm, the respective pressures are 510 mPa at 1 K and 
500 mPa at 300 K, giving a 2% reduction. 

3. When a = 500 nm, the pressures are 16.3 mPa and 15.2 mPa, giving 
a 6.7% reduction. 

4. For large gap widths the percentage corrections become much higher, 
though the pressures themselves are much weaker thus making experimental 
work more difficult. When a = 1 /im the respective pressures are 1.12 mPa 
and 0.96 mPa, giving a 13.9% reduction. 

It is of interest to compare the above findings with figure El in JU] • That 
figure shows, in the case of Au-Au, how the Casimir pressure varies with 
aT when the media are assumed nondispersive. The width is assumed to be 
a = 1 fim. The case temperature T = 300 K corresponds to aT = 0.131. One 
sees that in this case the agreement with our result above, T = 0.96 mPa, is 
reasonably good, if we put e = 3000. The broken line in figure [3] in JTU] gives 
the result when the modified ideal metal model is used in the calculation. As 
already mentioned, this is an idealized model, which assumes unit reflection 
coefficients for all but the TE zero mode: 

A = l, 50 = 0, A ro = S m =lform>l. (16) 

The MIM model corresponds to e — oo. Putting aT = 0.131 in the men- 
tioned figure El we see that T ~ 1.1 mPa. There is thus more than a 10% 
overprediction of the Casimir pressure following from the MIM model, at 
a = 1/iHi and T = 300 K, in comparison with our result 0.96 mPa above. 
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We recall again that at T = there is no distinction between a MIM 
model and an "ideal metal" model (IM), for which A m = B m = 1 for all 
m > 0. In the mentioned figure setting aT = 0, we obtain for an ideal 
metal T ~ 1.3 mPa at a = 1 /xm, which is considerably larger than the value 
1.12 mPa calculated above for T = 0. 

4 A temperature correction of experimental 
interest 

The calculated temperature correction at T = 300 K as compared with the 
T = case, although of fundamental interest, will be difficult to measure 
in practice. For practical purposes it can thus be better to focus on tem- 
peratures that are more realistic in the laboratory. In the following, as an 
example, we calculate the difference between the Casimir pressures at two 
temperatures T\ = 300 K and T 2 = 350 K, and let from now on AJF mean 
the pressure difference: 

A^= l^sool - |^350|. (17) 

This idea of testing the Casimir force seems to go back to Chen et al. |2*E] . 
and was elaborated upon also in [T2"] . 

Fi gures ITU1 and ITT1 show how the quantity (JT7j) varies with a in the interval 
a e [50, 1400] nm, for the same combinations of materials as before. Again, 
the differences between the materials are seen to be small. Taking Au-Au as 
an example, we see that AJ 7 = 2.0 mPa when a = 200 nm. It would perhaps 
be possible to measure a quantity like this. The experimental advantage we 
see of this kind of experiment is that only force differences are involved, for 
a given value of the gap width. Then there will be no need of measuring the 
absolute Casimir pressure itself, to an extreme accuracy. (Thermal expansion 
effects, of course, will have to be taken into account.) 
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Figure 10: Difference between Casimijgpressures at T\ = 300 K and T2 = 350 
K, equation (fTTjl . versus a for a G [50, 200] nm. 



Finally, figure shows how the relative magnitude of the difference 
Casimir pressure, AjF/jF 300 , for the same two temperatures. We see, in 
accordance with the behaviour above, that the relative temperature correc- 
tion is greatest when a is large. When a = 1.75 /im, the correction takes 
its maximum value, about 4%. Again, the experimental problem at large 
distances is that the forces themselves are so small. 
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Figure 12: Relative change of the Casimir pressure, between T\ = 300 K and 
T 2 = 350 K. 



5 Summary 

For similar and dissimilar plates, including the metals gold, copper, and 
aluminium, we have made accurate calculations of the Casimir pressure and 
have shown the results graphically. Basic elements in our calculations are 
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the Lifshitz formula; cf. equations (J3J) and (|7ji. together with realistic room- 
temperature values of the permittivities e(i(). For frequencies £ < 1.5 x 10 15 
rad/s; cf. equations (JTJ) and (J2J), the Drude dispersion relation is followed 
with great accuracy. 

We show the results at three chosen temperatures: (i) at T = 1 K, 
representing the T = case with good accuracy; (ii) at T = 300 K; and (hi) 
at T = 350 K. The low-temperature case is calculated numerically, without 
involving the modified ideal metal (MIM) model [TUj . 

It turns out that the differences between the Casimir pressures for the 
metals investigated here are small. From figures 6 and 10, for instance, it is 
seen that it is the case of Al-Al surfaces that gives the strongest forces. 

The most promising option for measuring the Casimir temperature cor- 
rection in practice seems to be to measure the pressure difference between 
two practically accessible temperatures in the laboratory. As figure 12 shows, 
the relative change of the Casimir pressure between the temperatures 300 K 
and 350 K are about 4% when a = 1.75 /im. The practical problem here 
is that the forces themselves are so small. For lower values of a the forces 
increase in magnitude, but the relative temperature corrections then become 
smaller. 
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